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HFAG Charm Mixing Averages 

B. Petersen 

CERN, CH-1211 Genve23, Switzerland 



Recently the first evidence for charm mixing has been reported by several experiments. To provide averages of 
these mixing results and other charm results, a new subgroup of the Heavy Flavor Averaging Group has been 
formed. We here report on the method and results of averaging the charm mixing results. 
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1. Introduction 

Almost since the discovery of charm mesons, mix- 
ing of D° — D mesons have been sought in analogy 
to the well-known K° — K mixing. Due to very effec- 
tive GIM suppression, the expected mixing rate in the 
charm system is much smaller than for kaons. Only 
very recently, the BaBar [l[ and Belle @ collabora- 
tions have reported the first evidence of charm mix- 
ing 1 . These results have renewed the interest from the 
theory community as the observed mixing rate could 
be caused by physics beyond the standard model or at 
least provide additional constraints on new physics. 

None of the mixing measurement have a significance 
above four standard deviations, but several have simi- 
lar precision for the mixing parameters. By combining 
the measurements we therefore obtain more precise 
values for the mixing parameters and exclude the no- 
mixing hypothesis with larger confidence. Combining 
the different mixing measurements is not completely 
straightforward, since not all measurements are sensi- 
tive to the same charm mixing parameters. 

The Heavy Flavor Averaging Group (HFAG) in 
2006 created a subgroup with the responsibility of pro- 
viding averages of charm physics measurements. One 
of the high priority tasks of this group is to combine 
the charm mixing measurements into world-average 
values for the fundamental mixing parameters. The 
first average assuming CP conservation was shown at 
FPCP \$. Besides those results, we here report the 
first results of combining mixing measurements where 
we allow for CP violation. 



2. Averaging Method 

Mixing is present in the D° — D system if the 
mass eigenstates, \D±) and I-D2}, differ from the flavor 

eigenstates, \D°) and \D ). Generally one can write 



|fi,2) = p\D°) ± q\D }. The variables of fundamental 
interest are the mass difference, AM = Mi — M 2 and 
decay width difference, Ar = Ti — r 2 between the 
two mass eigenstates. Traditionally, in charm mix- 
ing one uses the dimensionless variables, x — AM/T 
and y = Ar/2r, where T is the average decay width. 
CP violation in mixing or in the interference between 
mixing and decay would manifest itself as \q/p\ =/= 1 
and <fi = arg(q/p) ^ 0, respectively 2 . In addition CP 
violation could show up in the decay itself giving rise 
to decay mode dependent parameters. 

Most measurements do not directly measure (x,y). 
For instance in mixing measurements using D° — > 
K + ir~ decays there is a unknown strong phase, 5kk, 
so the results obtained are for x' = x cos 8k-k + 
ysin^io and y' = —x sin Sxir + V cos Sr-k- In the av- 
eraging procedure, we first combine measurements of 
the same parameters to obtain the more precise ob- 
servables. Most measurements are performed using 
likelihood fits and the combination is therefore per- 
formed by multiplying likelihood functions from each 
measurement and finding the new maximum. By com- 
bining likelihoods, correlations between observablcs 
and possible non-Gaussian tails are taken into ac- 
count. For measurements which are not using like- 
lihoods, we construct a likelihood using symmetrized, 
Gaussian uncertainties. To combine different types 
of measurements, the different combined likelihoods 
are recalculated as a function of (x, y, 5kv) minimiz- 
ing over any other variables. Sk-k is included since 
there is both a direct measurement 0] and by com- 
bining the D° — » K + ir~ measurement with the other 
measurement of x and y, one can also get a precise 
measurement of 5k When plotting confidence con- 
tours for (x, y) we minimize the likelihood over 5ktt- 

The combining of likelihood functions is currently 
only done for the CP conserving case. In principle it 
can be done also for the CP violating case by sim- 
ply having two more variables, \q/p\ and <j), in the 
final likelihood function. Unfortunately not all like- 
lihoods are currently available for the measurements 
which allow for CP violation. A simple combination 



1 Shortly after the CHARM2007 workshop additional results 
with evidence for charm mixing has been reported by the BaBar 
and CDF collaborations. In these proceedings we will summa- 
rize the status at the time of the workshop. 



2 The phase 
of decay mode. 



is for the moment assumed to be independent 
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is therefore performed by forming a x 2 of all measure- 
ments expressed in terms of the fundamental mixing 
parameters. The \ 2 assumes Gaussian errors, but cor- 
relations between observables in each individual mea- 
surement is taken account by using the full covariance 
matrix for each result. 



3. CP Conserving Averages 

The following averages were performed by adding 
log likelihoods from fits where CP conservation was 
assumed. 



3.1. Lifetime Ratio Average 

One can observe charm mixing by finding a differ- 
ence in the lifetime measured in decays to CP eigen 
states such as D° — > K + K~ and D° — > tt + tt~ and 
the mixed-CP decay D° — > tt + K~. We combine six 
results [1, H, [E 0> @> from such analyzes. All of 
these measure yep = tkixItuk — 1. In the limit of CP 
conservation one has yep — y. The average of the six 
measurements is yep = (1.12 ± 0.32) x 1CP 2 . This is 
3.5cr from the no-mixing hypothesis. As can be seen 
from Figure [TJ this average is mainly driven by the 
recent Belle measurement. 
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Figure 1: Measured j/cp values and the HFAG average. 



3.2. Mixing Rate Average 

Wrong-signed semileptonic decays provide a clean 
way of searching for charm mixing, but the mea- 



surements are only sensitive to the integrated mix- 
ing rate Rm = {x 2 + y 2 )/2. Four measurements 
[To . [TTI . [T2I Il3j are combined and give an average 
of R M = (1.7 ± 3.9) x 1CT 4 . In addition to the 
semileptonic decays, Rm can also measured in the 
analysis of fully hadronic decays. The semileptonic re- 
sult is therefore combined from two hadronic analyzes 
(l4l HU and in addition an analysis of tagged decays at 
the ^(3770) Q. The combination is illustrated in Fig- 
ure[2]and gives an average value of Rm = (2.1±1.1)~ 4 . 
In the transformation to a likelihood in (x,y), we ig- 
nore the non-physical region of Rm < 0. 
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Figure 2: The mixing rate from measurements using 
semileptonic D° decays are averaged with results from 
multi-body hadronic charm decays. 



3.3. (x, y) Average 



One can measure x and y directly using a time- 
dependent Dalitz plot analysis of D° — > K*gir + ir~ de- 
cays. Two measurements [H, [I3| have been published 
and these have been averaged by HFAG and gives 
x = (8.1 ± 3.3) x 1CT 3 and y = (3.1 ± 2.8) x 1CT 3 . 
Combining this average with the averages above for 
Rm and ?/cp using likelihoods mapped as a func- 
tion of (x, y) we obtain x = (9.2 ± 3.4) x 1CT 3 and 
y = (7.0 ± 2.2) x 10~ 3 . Contours of the combined 
likelihood function at the levels corresponding to 1 
to ha confidence levels are shown in Figure [3] Note 
that the confidence levels shown correspond to two- 
dimensional coverage probabilities of 68.27%, 95.45%, 
etc., and therefore 2AlnL = 2.30,6.18, etc. 
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Figure 3: Confidence level contours in the mixing pa- 
rameters (x,y) from the combination of yep, Rm and 
(x,y) from the time-dependent Dalitz analysis of D° — > 



3.4. Averages for D° — ► K + tv~ Decays 

As mentioned above, one can measure x' and y' 
using the doubly-Cabibbo suppressed (DCS) decay 
D° — > AT + 7r — . The likelihood functions are available 
for two measurements [H [HI of this type. These are 
combined and gives the averages x' 2 = (—0.1 ± 2.0) x 
10~ 4 and y' = (5.5_tf;f) x 10~ 3 . The corresponding 
likelihood contours are shown in Figure [4j 
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Figure 4: Confidence level contours from the combination 
of BaBar and Belle measurements using D° — * K + tv~. 
The wiggles in the 4 and 5a contours for x' 2 < is a 
binning effect. 



other mixing results in Section 13.31 which do not de- 
pend on Sk-tt- An additional constraint comes from 
a CLEO-c measurement [4[ of cosSktc = 1-09 ± 0.66, 
where a small dependence on x and y is ignored in 
the combination. Figure [5] shows the likelihood con- 
tours in (x,y) after minimizing over 5ktt- The region 
around the central value is almost unchanged with re- 
spect to the result without the D° — > K + n~ decays 
(Figure[3]). This is also reflected in the over all average 
for x and y which are 

x = (8.7t 3 ;°) x 10- 3 , 
y = (6.6 ±2.1) x 10~ 3 . 

The D° — » K + tt~ measurements do not contribute 
much to the central value, because of the poorly 
known phase 5ktt- However they do help exclude the 
no-mixing hypothesis and cause the dip seen in the 
contours close to (x,y) — (0,0). At (x,y) — (0,0) we 
obtain 2A In L = 37 with respect to the minimum. 
This corresponds to a significance of the combined 
mixing signal of 5.7o\ 
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Figure 5: Confidence level contours from combining all 
mixing measurements under the assumption of CP con- 
servation. 

The combination also gives an improved value for 
Sktt- This can be seen in the projection of the like- 
lihood after minimizing over x and y in Figure [HI 
The combination gives Sktt(= 0.33lg;2g)rad. With- 
out the CLEO-c measurement of Sktt, there would be 
an equally good second minimum at Sktt = — 2.17rad. 



3.5. World Average 

The combined likelihood for D° — > K + ir~ decays 
can be expressed as a function of (x, y, 8k-k) ignoring 
the part with x' 2 < 0. This likelihood can be com- 
bined with the likelihood from the combination of the 



4. CP Violating Averages 

Measurements of charm mixing can be done with- 
out assuming CP conservation by fitting D° and D 
mesons as separate samples. Most of the measure- 
ments above have done that and we therefore can com- 
bine those to also provide constraints on the CP vi- 
olating parameters. When allowing for CP violation, 
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Figure 6: Log-likelihood function for the strong phase, 
8k tv, from combining all mixing measurements. 



The mixing parameters are almost unchanged with 
respect to the CP conserving average. This is also seen 
from the confidence levels shown in Figure [7] One can 
also draw the 1 to 5a confidence level contour for <fi 
versus \q/p\ using A% 2 . This is shown in Figure[H The 
no-CP violation hypothesis is seen to lie well within 
the la contour. 

The combined results for the D° — > K + ir~ param- 
eters are 

5 K * = Q.33±°;|>rad, 

R D = (3.35 ±0.11) x lfr 3 , 

A D = (-0.8 ±3.1). 

There is little change in Sk* with respect to the CP 
conserving average and no evidence for direct CP vi- 
olation as Ad is consistent with zero. 



the measured parameters are related slightly differ- 
ently to the mixing parameters. For the lifetime ratio 
measurements, one has 

2?/cp = (\q/p\ + \p /q\)y cos (f>- (\q/p\ - \p/q\)xsin(p, 
2A T = (\q/p\ - \p/q\)y cos <j>- (\q/p\ + \p/q\)x sin (j), 



where At is the measured relative lifetime difference 
for D° -> h+h~ and if -> h+h~. For D° K + it~ 
decays, the x' and y' measured for D° and D are 
related as follows 



,.'± 



,.'± 



l±A 



M 



1 T A M 
1 ± Am 



1/4 



(x' cos 4>±y sin < 



l T A 



M 



1/4 



{y' cos 4>-fx' sin < 




where A M = ■ For D° -> K%^- de- 

cays the measurement directly gives x, y, \q/p\ and 
cj), while for the Rm analysis the results are not sepa- 
rated and therefore just measure Rm = {x 2 + y 2 )/2. 
The measurement of 5ktt from CLEO-c is not done 
separately for D and D mesons and is not included 
in the combined result allowing for CP violation. 

In total 22 measurements are combined in a x 2 -fit 
to extract seven parameters, the four mixing and CP 
violation parameters, x, y, \q/p\ and <p, and three char- 
acterizing D° —> K + ir~, namely Skit, the DCS rate 
Rd , and the direct decay rate asymmetry Ar> ■ The fit 
gives x 2 = 14-4 and the following mixing parameters 

x = (8A±H) x lO" 3 , 

y = (6.9 ±2.1) x 10~ 3 , 

\q/p\ = 0.88l°; 23 , 

= (-0.09l8:g) rad. 



Figure 7: Confidence level contours for (x, y) from com- 
bining mixing measurements with CP violation allowed. 
The dashed blue curve shows the lcr contour from the CP 
conserving case for comparison. 



5. Summary 

Evidence of charm mixing has been reported from 
several experiments in the last year. A new subgroup 
of HFAG has performed an average of these and other 
existing charm mixing results. The combined result 
has a signal significance in excess of 5 standard devi- 
ations and gives the mixing parameters 

x = (8.4l 3 ; 2 ) x 10- 3 , 
y = (6.9 ±2.1) x 10~ 3 . 
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Figure 8: Confidence level contours for (\q/p\,(f>) from 
combining mixing measurements with CP violation al- 
lowed. 



CP violation parameters have also been combined and 
gives 



\q/p\ = 0. 



d+0.23 
'-0.20' 



= (-0.09t°;i£)rad. 

This is fully consistent with no CP violation being 
present in charm mixing. HFAG intends to period- 
ically update these averages as new results become 
available in order to provide the most precise mixing 
parameters to the community. 



References 

[1] B. Aubert 
tion], Phys. 



et al. [BABAR Collabora- 
Rev. Lett. 98, 211802 (2007) 



arXiv:hep-ex/0703020 



[3 
[4 

[5; 

[G 
[7 

[8 

[9 
[10 

[11 

[12 
[13 

[14 

[15 
[16 

[i7: 

[18 



A. J. Schwartz, In the Proceedings of 5th Fla- 
vor Physics and CP Violation Conference (FPCP 
200 7), Bled, Slov enia, 12-16 May 2007, pp 024 
|arXiv:0708.4225l [hep-ex]] . 

D. M. Asner et al. [CLEO Collaboration], 
Int. J. Mod. Phy s. A 21, 5456 (2006) 
|arXiv:hep-ex/0607078| . 

E. M. Aitala et al. [E79 1 Collaboration], Phys . 
Rev. Lett. 83, 32 (1999) |arXiv:hep-ex/9903012| . 
J. M. Link et al. [FOC US Collaboration!, Php . 
Lett. B 485, 62 (2000) [arXiv:hep-ex/0004034l . 
S. E. Csorna et al. [CLEO Collabora- 
tion], Phys. Rev. D 65, 092001 (2002) 
|arXiv:hep-ex/0111024| . 

B. Aubert et al. [BABAR Collabora- 
tion], Phys. Rev. Le tt. 91, 121801 (2003) 
[arXiv:hep-ex/0306003| . 

K. Abe et al. [Belle Collaboration], Phys. Rev. 
Lett. 88, 162001 (2002) [arXiv:hep-ex/0111026l . 
E. M. Aitala et al. [E791 Collabora- 
tion], Phys. Rev. L ett. 77, 2384 (1996) 
arXiv:hep-ex/9606016| . 

C. Cawlfield et al. [CLEO Collabora- 
tion], Phys. Rev. D 71, 077101 (2005) 
|arXiv:hep-ex/0502012| . 

U. Bitenc et al. [Belle Collaboration], Phys. Rev. 
D 72, 071101 (2005) [arXiv:hep-ex/0507020] . 
B. Aubert et al. [BABAR Collaboration], Phys. 
Rev. D 76, 014018 (2007) [ arXiv:0705.0704l [hep- 
ex]]. 

B. Aubert et al. [BABAR Collabora- 
tion], Phys. Rev. Le tt. 97, 221803 (2006) 
|arXiv:hep-ex/0608006l . 

B. Aubert et al~ [BABAR Collaboration], 
|arXiv:hep-ex/0607090| 

D. M. Asner et al. [CLEO Collabora- 
tion], Phys. Rev. D 72, 012001 (2005) 
[arXiv:hep-ex/0503045| . 

L. M. Zhang et al. [Belle Collaboration], Phys. 
Rev. Lett. 99, 131803 (2007) |arXiv:0704.i7j00] . 
L. M. Zhang et al. [BELLE Collabora- 
tion], Phys. Rev. Le tt. 96, 151801 (2006) 
|arXiv:hep-ex70601029| . 



[2] M. Staric et al. [Belle Collaboration], Phys. Rev. 
Lett. 98, 211803 (2007) [arXiv:hep-ex70703 036 . 



